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1.  Introduction 


Understanding  and  using  Double  Negative  (DN)  materials — negative  a{ co)  and  /j(co)  for  some 
frequencies,  in  the  design  of  devices  using  electromagnetic  waves  remains  a  challenge.  This 
topic  has  become  of  great  interest  since  these  man-made  materials  became  popular  this  decade. 
Many  of  us  had  difficulties  understanding  propagation  because  we  had  not  experienced  DN 
materials  before  and  we  had  to  except  the  idea  of  a  negative  phase  velocity  as  a  consequence  of 
our  equations.  The  concept  of  a  positive  phase  velocity  appears  to  have  always  been  sacrosanct. 
It’s  interesting  to  note,  however,  that  scientists  could  accept  the  idea  of  a  phase  velocity  going 
faster  then  the  speed  of  light  because  invariably  a  real  physical  situation  associated  with  this 
occurrence  could  be  found. 

Here  we  derive  the  propagation  equations  in  DN  materials  in  a  way  that  differs  from  previous 
derivations.  We  propose  that,  for  any  real  material,  the  energy  of  the  wave  must  decrease  as  it 
moves  from  its  source  of  energy.  Energy  losses  per  unit  length  of  travel  are  finite,  and  can  be 
extremely  small.  Propagation  in  a  lossless  DN  media  is  found  as  the  mathematical  limiting 
solution  of  an  extremely  small  energy  loss  per  unit  length  for  real  s(co)  and  ju(co).  We  end  up 
with  the  same  conclusions  as  the  other  authors  using  nothing  more  than  an  algebraic  solution  to 
Maxwell’s  equations  without  the  need  to  even  mention  phase  velocity;  negative  phase  velocity  is 
a  consequence  of  our  solution. 

Why  was  it  necessary  to  develop  this  derivation  when  the  result  was  already  available?  We’ve 
gone  through  the  trouble  because  we  plan  to  examine  complex  multi-dimensional  structures  and 
do  not  want  to  confront  potential  ambiguities  about  which  signs  to  choose  in  taking  square  roots, 
etc.  Having  a  guiding  principle  to  go  by  is  a  big  help.  In  this  report  we  apply  the  basic 
principles  to  reflection  and  refraction  of  an  electromagnetic  waves  between  a  material  having 
positive  a{co)  and  ju(co),  and  one  having  negative  a{ co)  and  ju(a>). 

In  the  next  section  we  develop  the  propagations  using  the  principles  just  mentioned.  Use  of  the 
equations  to  solve  the  aforementioned  reflection/refraction  problem  is  rendered  in  section  3. 
Concluding  remarks  are  rendered  in  section  4  and  related  references  are  listed  in  section  5. 


2.  Electromagnetic  Propagation  Model  for  DN  Materials 


For  a  DN  material  Maxwell’s  equations  are 

V*B  =  0  (1) 
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V  •  D  =  0 
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Using  the  formulas 
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we  write  equations  5  and  6  in  component  fonn 
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The  constitutive  equations  are 


D{t)  =  e0E(t)  +  P(t) , 

B  =  MH(t)  +  ^(0)  =  Hflit)  +  B0M(t) , 


(3) 
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A  dynamic  relationship  is  assumed  between  the  polarization,  P(t) ,  and  E{t) ,  and  between  the 
magnetization,  M(t ) ,  and  Flit) .  An  example  of  dynamic  equations  for  the  polarization  and 
magnetization  is  the  set  provided  by  Smith  (16): 

d2P  „  dP 


dt2  +rE  dt  +0)eoP-£o<E 


d-M  „  dM  2  -  2  - 

~dF+ r*-dT+a>m0' - 


(9) 

(10) 


The  system  of  equations  is  solved  with  the  fonnula  (/dt)  =jco  and  using  the  transform  pair 

—  1  +co 

F(t )  =  —  }  zxp(j  cot)F  (co)dco 

In  -oo 

^  +oo  __ 

F(co)=  J  Qxp(-jaX)F(t)dt 

-oo 

The  function  F(t)  represents  either  P(t) ,  E(t) ,  M(t),  H(t) ,  or  J .  We  have 

P  =  £0ZeE 
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For  illustrative  purposes  we  use  the  following  fonn  for  %E(cd),  although  this  is  never  necessary  in 
the  calculations. 
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B(co)  =  ju0(  1  +  %(co))H(co)  =  ju0jur(co)H(co )  =  ju{co)H(co) 
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As  done  for  %F(co),  we  use  an  illustrative  expression  for  zn( co) 
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In  Fourier  transfonn  space  equations  5  and  6  now  become 
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Now  let  us  consider  propagation  in  the  positive  z-direction  where  there  are  no  current  sources. 
Let  J  =  0  for  z  >  0.  Thus  for  z  >  0  the  magnitudes  of  the  fields  must  decrease  as  z  increase  for  a 
lossy  medium  because  all  the  sources  of  energy  lie  in  the  region  z  <  0. 
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For  brevity  we  choose  a  TEM  wave  with  Ex  and  Hy .  We  also  assume  d/dx  =  0  and  d/dy  =  0. 


Equations  29  and  30  then  become 
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We  seek  solutions  of  the  form 

Ex  =  E  cxp(-  jkz) 

(35) 

Hv  =  H  exp(-jkz) 

(36) 

and  substitute  these  expressions  into  equations  3 1  and  32.  We  get  dldz  =  -jk  and 

-  jkE  =  -  jcOju(co)If 

(37) 

jkH  =  jcoa{(o)E 

(38) 

A  solution  occurs  only  when 

k 2  =  ora(oj)jU(oj) 

(39) 

All  of  the  equations  presented  so  far  have  been  derived  before  in  various  forms  in  most  books  in 
electromagnetic  wave  propagation.  Solutions  for  k  have  been  addressed  for  three  cases:  a  and  // 
are  both  positive,  a  is  positive  and  //  is  negative,  and  a  is  negative  and  //  is  positive.  We  are  now 
considering  the  case  where  both  a  and  //  are  both  negative.  We  show  that  the  solution  is  a 
straightforward  extension  of  the  three  previous  cases  to  this  new  regime,  and  needs  no  additional 
mathematical  or  physical  assumptions  other  than  using  the  fact  that  propagation  material  is  lossy, 
albeit  extremely  minute. 

Before  going  further,  a  few  comments  about  using  the  foregoing  set  of  equations  is  appropriate. 
First,  we  are  dealing  with  casual  functions,  whose  mathematical  properties  are  well  known. 

Being  casual  we  need  only  deal  with  the  co  >  0  part  of  the  complex  plane.  From  equations  20, 
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28,  33c,  and  34c  we  then  see  that  Pe{co)  and  Ph(co)  are  both  positive.  In  addition,  we  also  notice 
that  «/■(  co)  and  au(co)  are  even  functions  of  co,  namely  they  are  functions  of  of.  These  are 
important  properties  in  the  analysis. 

Also,  even  though  our  analysis  uses  specific  forms  for  ai(co),  C(h(co),  /f  ( co),  and  Puico),  the 
aforementioned  symmetries  of  the  functions  are  general.  They  apply  to  any  casual  function. 

Using  equations  33a  and  34a  in  equation  39  we  get 

k  —  ko  —  J Pe  )(^//  — ./ Ph )  —  ko  —  PEPH )  —  jk()  {ctEPH  +  ci n [3 F )  (40) 

2 

i  2  2  ^  /  a  a  \ 

k0  =  co  £0//0  =  y  (41) 

c 

where  c  =  (1  /  is  the  velocity  of  light  in  vacuum.  The  right  hand  side  of  equation  40  is  a 

complex  number  and  therefore  k  will  be  a  complex  number.  Let  us  write  k  as 

k  =  k0(rj  +  j 4)  (42) 

where 

7  +  j£  =  ±((oceoCh  ~  PePh  )  ~  MePh  +  ofA  ))'  2  (43) 

When  equation  42  is  used  in  equations  35  and  36  we  get 

Ex  =  E  exp {-jkz)  =  E  exp (-jk0?i  z)  exp(kng  z)  (44) 

H  =  H  exp  {-jkz)  =  H  cxp(-jk(]//  z)  cxp(  k(jc  z)  (45) 

If  the  fields  are  to  diminish  as  z  increases,  then  c,  must  be  negative.  This  is  how  the  +  or  -  sign  is 
chosen  in  equation  43. 

For  orientation  let  us  consider  the  traditional  case  first.  This  is  the  one  where  the  a-terms  are 
positive  and  are  much  larger  than  the  /?  terms.  Equation  43  is  then  approximated  by  the  equation 
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Since  the  (/fefi-terms  are  each  «1  we  can  use  the  approximation 
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(48) 
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If  the  “+”  is  used  for  the  right  hand  side  of  equation  48,  then  £  will  be  negative  and  satisfy  the 
conditions  of  equations  44  and  45.  This  is  therefore  the  correct  decision.  We  then  get 
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Because  cfe(ru),  cch(co),  Pe(co),  and  fn( co)  are  all  positive,  rj  '  is  positive  and  <f+)  is  negative. 
Equations  49  and  50  are  what  we’re  used  to  seeing.  As  we  shall  show,  this  result  leads  to  the 
notion  of  a  positive  phase  velocity  in  the  same  direction  of  energy  flow. 


When  equations  49  and  50  are  inserted  into  equations  44  and  45  we  immediately  notice  that  the 
phase  part,  exp(-/'A:o//z),  is 


exP (~jk0ri{+)  z)  =  exp(-  jk0(aEaH)V2  z)  =  exp 
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(51) 


Since  the  time  dependence  is  exp (jcot)  we  see  that  equation  5 1  defines  a  frequency  dependent 
phase  velocity, 
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The  attenuation  is  given  by  exp (k0<^z)  which  is 
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The  foregoing  system  of  equations  that  we  have  just  analyzed  is  called  a  right-hand  (RH)  system. 
In  our  construct  it  is  comprised  of  the  following  three  vectors: 
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foregoing  equation  ax ,  a  ,  and  az  are  the  unit  vectors  in  the  x-,  y-,  and  z-directions.  Since 
positive  the  three  vectors  in  equation  54  follow  the  RH  rule,  that  is 
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E  x  H  =  EHa 
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The  RH  rule  is  natural  for  cases  where  the  phase  velocity  points  in  the  positive  z-direction.  This 
is  also  physically  comforting  since  the  root  mean  square  Poynting  vector,  Srms  also  point  in  the 
positive  z-direction.  Using  equations  45  and  46  we  have 

Srms=^Re(ExH*)  =  ^:RQ(m*)exp(2k0?+)z)  (56) 


Using  equation  50  we  write 

Srms  =  ^-Re(£  x  H*)  =  X-az  Re(EH*)exp(-k0/+)z ) 
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The  last  remaining  step  to  calculate  Srms  is  to  compute  Re(£7/!)  from  either  equation  37  and  38. 
Equation  38  gives 


k 


(59) 


Using  equation  33a  for  s(a>)  and  equation  42  for  k (co)  we  get 
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In  the  foregoing  equation  Z0  =  ^ju0  /  e0  is  the  impedance  of  free  space, 


(aEaH)U2  -  —  jy{+) 

g(+)  = - - SP(+)exp(y^+)),  (61) 

-  JPE 

P{  1  and  1  are  the  amplitude  and  phase  of  g(+)  respectively,  and  Z(+)  =  g(+)Zo  is  the  wave 
impedance.  Since  all  the  parameters  in  g  are  positive  the  phase,  tP+\  is  less  than  n/2.  The  real 
part  of  g(+)  is  positive. 
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Thus, 


SW  =  a 


R  e{EH*)  = 

HI 


cos1/7 


(+) 


R(+,Zn 


cos1/7 


(+) 


2  R(+)Z, 


-exp(-k0y(+’z) 


(62) 

(63) 


As  we  see  from  equation  55  the  three  vectors:  E,  H  and  Snm  form  a  RH  system  since  Srms 
points  in  the  positive  z  direction. 


Now,  let  us  see  what  happens  at  the  other  extreme,  when  both  «£•  and  cch  are  negative.  This 
happens  for  the  set  of  radian  frequencies,  { (o\ ,  that  satisfies  the  conditions 


We  write 


«£(®)  =  1  +  j;(/y)<0 

(64) 

M®)  =  i  +  zU®)<° 

(65) 

| 

(66) 

<2^  =  —  |<2^  | 

(67) 

±(( \aE \aH  \~PePh)+  j(\aE  \PH  +  \aH  \Pe  ))‘/2 

(68) 

Let’s  first  evaluate  equation  68  when  the  /Tterms  are  much  smaller  than  the  a-terms.  This  is 
what  we  did  before  and  is  the  usual  case  of  interest.  We  now  get 


V  +  J  Z  =  ± 


(«/ 


|l/2  | 


a 


1 1  /  2 


H 


)  +  i 


\ac 


1 1/2 


|l/2 


Pt!  |  g// 


1 1/2 


|l/2 


yy 


(69) 


Again,  imposing  the  condition  that  c  be  negative  to  ensure  energy  decay  in  the  positive  z,  we 
now  require  that  we  take  the  negative  sign  of  the  right  hand  sign  of  equation  69.  In  lieu  of 
equations  49  and  50  we  now  have 

7(-}  =  -(k|k|)1/2  (VO) 


k}=-- 

r\ 


r  1V2  n  r  11/2  a 
[ aE  Ph  ,  [aH  Pi 


+  ■ 


M"  hJ 


V  L  H 


/  2 

e  j  y 


-y 


(-) 


z*  * = 


[a,-]  7/,  ,  [«„]'  7/' 


kb  kf 


V  L^// 


/2 

£j  y 


(71a) 

(71b) 
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Notice  that  £  )  and  df+)  are  both  negative:  compare  equation  71  with  58. 

When  equations  70  and  71  are  inserted  into  equations  44  and  45  we  immediately  notice  that  the 
phase  part,  exp(-jk0r/z),  is 


exP(-7^o77<_)  *)  =  exp(y^0(|«£|«//|)1/2  z) 


=  exp 


f  .  \ 

v(-) 

V  vp/!  7 


(72) 


Since  the  time  dependence  is  Qxpijcot)  we  see  that  equation  72  defines  a  negative  frequency 
dependent  phase  velocity, 


VpH  =  ~ 


(Mk|)1/2 


(73) 


The  attenuation  is  given  by 


exp(£0£(  1 z)  =  exp 


k0z 


(\  1 1/2  „  I  ,1/2  A 

l^fel  A/  +  |a//|  Pe 


V  1 


,1/2 


ac 


,1/2 


exp 


k0Z  y(~) 


(74) 


As  before  we  have 


~  _  O)S(O))  ~  _ 


«£(«) 


k  k0(fi^  + 

but  because  of  equation  64  this  time  we  need  to  use  the  expression 

£(a)  =  £fl(a/:(oj)  -  jpE(oj))  =  -e0(\aE\  +  jpE) 


(75) 


(76) 


From  equations  70  and  71  we  have 


//->  =-  (a£  a„)1,z  +^r 


.  Z„<-> 


7 


(77) 


Inserting  equations  76  and  77  into  equation  75  and  using  previous  definitions  gives 

B=-* 


Z(^ 

z(-)=g(-)z0 

1 


g(~’  = 


"*■  j Pe 


=  P{  )  exp (y  IP1  () 


(-)> 


^rms 


E  2  COS^*  ’ 


2Z0,P 


(-) 


-expi-ky^z) 


(78) 

(79) 

(80) 

(81) 
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Again  we  see  from  equation  8 1  that  the  Poynting  vector  points  in  the  positive  z-direction,  and 
therefore,  Srms ,  E  ,  and  H  form  a  right  handed  system.  The  fact  that  this  situation  is 

accompanied  with  a  negative  phase  does  not  violate  any  law.  Both  cases,  the  double  positive 
and  double  negative  a e  and  a h  are  described  within  the  same  mathematical  framework  for  a 
plane  wave  traveling  in  any  direction. 

Let  Qe  be  the  direction  of  the  electric  field  and  Ou  be  the  direction  of  the  magnetic  field.  We 
have 


Qe  •  Qh  =  0 

(82) 

Qs  =  Qe  x  Qh 

(83) 

Qs  •  Qe  =  Qs  •  Qh  =  0 

(84) 

E  =  E0Qe  cxp(-  /A0//Q  v  •  r)  exp {-ka]Els  •  r  ) 

(85) 

~  E  ~ 

H  =  -y- El H  exp(-/l0?;Qs.  •  r)  exp(-k0^Q5  •  r) 

(86) 

For  double  positive  materials  use 

7(+)  =  (aEaH)U2  ,  /(+)  = 
and  for  double  negative  materials  use 

7(_)  =-(K||aj?|)1/2»  = 


rvU2  nvl 

V  UH  UE  J 


_  y(+) 


z  =  z 


(87) 


k]I/2^  +  kf2A 


|l/2  o  \ 
E 


Kt2  kf2 


-  y(-) 


Z  =  Z 


(88) 


3.  Reflection  and  Refraction  Involving  a  DN  Material 


In  this  section  we  demonstrate  the  unique  features  of  reflection  and  refraction  involving  a  DN 
material.  In  the  three  cases  considered  in  figures  1  through  3  the  electric  field  is  in  the  plane  of 
incidence.  The  behavior  of  figure  1  is  well  known;  we  include  this  for  the  reader’s  orientation. 
We  demonstrate  the  results  shown  in  figures  2  and  3  from  the  solution  of  the  equations. 
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Figure  2.  Reflection  in  air  off  a  negative  sand  negative  u  material. 
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Figure  3.  Reflection  in  a  negative  s  and  negative  u  material  off  an  air  interface. 


In  all  cases  considered  in  figures  1  through  3  equations  29  and  30  apply.  They  are  simplified 
according  to  equation  89. 


r 


V 


rdE,  dEv A 


By 

8z 

J 

{  BEX 

8EZ 

l  & 

8x 

r  be  v 

8EX 

\ 

V 

8x 

8y 

J 

8HZ  8Hy" 

8y  8z 
r  8H ,  8H_ 


=  -jcoju(co)Hy 

=  - j  coju(co)H  _  =  0 
8HV 

— -  =  ja>£(co)L 
8z 


v  8z  8x  j 


8H  8HX 


=  jco£(co)Ev  =  0 


8x  8y 


8H 


8x 


-  =  j0)£((O)Ez 


dH'=jcoeUo)dS 


8z 


8z 


X 


(89a) 

(89b) 

(89c) 

(89d) 

(89e) 

(89f) 

(90a) 
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d  HV  ■  < 

—^  =  JCDS(cd)—^ 
ox~  ox 

Inserting  equations  90a  and  90b  into  equation  89b  gives 

dHv  dHv 

— {-  +  —t  =  -cdju(cd)s(cd)H 
oz  ox 

For  the  case  where  ju(co)  and  s( co)  we  get  the  familiar  results 


(90b) 


(90c) 


A  exp  [-y  ( kxx  +  /czz)]  =  Aexp-  —  (Q»r) 

(91) 

CDS 

(92) 

l 

CDS 

(93) 

r  =  xax  +  zaz 

(94) 

Q  =  sin  cpcix  +  cos  cpaz 

(95) 

1 

V  =  , - 

(96) 

In  the  foregoing  equations  ax  and  a,  are  unit  vectors  in  the  x-  and  z-  directions  respectively,  Q  is 
the  direction  of  propagation,  sin  (p  =  Q»ax,  and  cos  cp  =  Q  •  a. .  For  a  wave  traveling  in  the 
positive  x-  and  positive  z-direction  we  have 


,  co  . 
kx  =  — sin  cp 

V 

(97) 

CD 

k,  =  — cos  cp 

V 

(98) 

—  H  =  Z  cos  cpH 

CDS 

(99) 

— H  =-Z  sin  (pH 

CDS 

(100) 

N 

II 

(101) 

Solution  for  reflection  in  air  off  a  positive  e  and  positive  p  material 
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From  figure  1  we  have  the  following 

Q,  •  r  =  x  sin  cp1  +  z  cos  cp, 
ClR  •r  =xsin  <pR  -  z  cos  cpR 
Qr  •  r  =  x  sin  cpT  +  z  cos  cpT 


Using 


~  1  dHv 

E_=  y 


E..  = 


jcos  Gx 
jcos  Gz 


H  j  =  El  exp-  —  j(Qj  •?)  =  H  exp-  —  j(x  sin  cpI  +  z  cos  cpt ) 


co 


ExI  =  Z0  cos  cp,E[  exp - j(x  sin  (p,+z  cos  cp, ) 


CO 


E  j  =  -Z0  sin  (ptE[  exp - j(x  sin  cpI  +  z  cos  (p, ) 


Hv,r  =  A  exp-  —  j(ClR  »r)  =  A  exp-  —  j(x  sin  <pR-z  cos  <pR  ) 


CO 


ExR  =  -Z0  cos  cpRA  exp - j(x  sin  cpR-z  cos  cpR  ) 


co 


E  R  =-Z0sin^Zexp - j(xsmcpR  -zcos(pR) 


HvT  =B exp-  —  j(ElT  •r)  =  B exp-  —  j(x sin cpT+z cos cpT ) 

V  V 


CO 


Ex  T  =  Zm  cos  cpTB  exp - j(x  sin  cpT  +  z  cos  cpT ) 


Ev  T  =  -Zm  sin  <pTB  exp-  —  j{x  sin  cpT  +  z  cos  cpT ) 

v.„ 


1 


(a) 


1 


v_  = 


V 


(b) 


(102a) 

(102b) 

(102c) 

(103a) 

(103b) 

(104a) 

(104b) 

(104c) 

(105a) 

(105b) 

(105c) 

(106a) 

(106b) 

(106c) 

(107) 


15 


Z0  =  J —  (a)  Zm=J—  (b) 


(108) 


In  the  foregoing  expression,  H  is  the  magnitude  of  the  incident  magnetic  field,  A  is  the 
magnitude  of  the  reflection  coefficient  and  B  is  the  magnitude  of  the  transmitted  wave.  These 
quantities  are  determined  by  matching  boundary  conditions  at  z  =  0 


Hv  ,(x,z  =  0)  +  Hv  R(x,z  =  0)  =  Hv  T(x,z  =  0) 


(109a) 


Hex p  y(*sin^7)  +  Aex p  -  —  j(xsirupR)  =  B ex p  ~  —  j(xsin(pT)  (109b) 


ExI{x,z  =  0)  +  ExR{x,z  =  0)  =  E  T(x,z  =  0) 


(110a) 


Z0  cos (pjH  exp  j{xsm(pI)  -Z0cos^Zexp  j(xsm(pR)  =  Zm  cos (pTB  exp  ~  —  j(xsin(pT) 

L  v0  J  L  V0  J  L  V,n 

(110b) 

The  solution  of  equations  109  and  1 10  is  available  in  standard  electromagnetic  texts.  We  have 


(Pr=(P, 


sm  (Pi  =  sin  (pT 


(111a) 

(111b) 


H  +  A  =  B 


Z0  cos  <pj  (H  -  A)  =  Zm  cos  <pTB  =  Zm  cos  (pT  ( H  +  A) 


(111c) 

(11  Id) 


A_Z0  cos  (pI-Zm  cos  cpT  H 
Z0  cos  cpl  +Zm  cos  (pT 


(1 1 le) 


B  = _ 2Z0  cos  cp, _ jj. 

Z0  cos  cp,  +  Zm  cos  cpT 


(Hlf) 


Solution  for  reflection  in  air  off  a  negative  8  and  negative  ju  material 

The  solution  for  this  case  is  readily  found  using  the  discussion  for  the  double  case  combined 
with  the  basic  conclusions  deduced  in  section  2.  The  basic  behavior  for  propagation  for  a  DN 
material  is  summarized  from  equations  84  to  88  applied  to  the  loss-free  model  considered  here. 
Restating  them  in  the  notation  of  this  section  we  have 


f2s  •  Qe  =  E2s»  Z2h  =  0 


(112) 


E  =  E0ClEexp(j  —  ns*r) 


(113) 
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(114) 


~  —  ft)  — 

H  =  -=-ClH  exp(j  —  Q5  •  r) 

Z  v 

m  m 

Equations  113  and  114  have  a  negative  phase  velocity. 

Suppose,  for  example,  that  we  assume  the  transmitted  wave  looks  like  the  one  for  the  double 
positive  case  just  considered — that  is,  figure  1  is  valid.  What  would  the  solution  for  the  system 
of  equations  look  like?  The  behavior  for  the  incident  (equation  104)  and  reflected  terms 
(equation  105)  remains  unchanged.  However,  the  equation  for  the  transmitted  terms  is  now 


co 


H  T  -B  exp — /(xsin  cpT  +  z  cos  <pT) 


co 


Ex  T  =  Zm  cos  (prB  exp  —  j(x  sin  cpT+z  cos  cpT ) 


CO 


Ez  T  =  -Zm  sin  cpTB  exp  —  j{x  sin  cpT+  z  cos  cpT ) 


(115a) 

(115b) 

(115c) 


At  the  boundary,  z  =  0,  the  spatial  variation  of  equation  115  goes  like 

Hv  T  — >  exp  —  j(x  sin  cpT  ) 
v 

m 

co 


"x,T 


exp  —  y(*sin  (Pj) 


(116a) 

(116b) 


By  comparing  the  foregoing  behavior  with  incident  and  reflected  behavior  at  z  =  0,  which  go  as 


exp 


co 


- y'(*sin^) 


,  we  see  that  there  is  no  way  a  transmitted  wave  traveling  in  the  direction 
shown  in  figure  1  can  satisfy  the  boundary  conditions. 

Now  let’s  consider  the  transmitted  wave  shown  in  figure  2.  For  this  case  we  have 

ClT  •  r  =  -x  sin  cpT  +  z  cos  cpT 

Again,  using  the  concept  of  negative  phase  velocity  we  now  have 


(117) 


H  T  =  B  exp  —  j(-x  sin  cpT  +  z  cos  (pT ) 

V 

m 

(118a) 

=  Zm  cos  cpTB  exp  —  j(-x  sin  cpT  +  z  cos  (pT ) 
v 

m 

(118b) 

=  Zm  sin  cpTB  exp  —  j(-x  sin  cpT  +  z  cos  cpT ) 

Vm 

(118c) 
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Repeating  the  same  procedure  as  in  equations  109  and  1 10  we  have 


//exp 

~—j(x  sirup,) 

+  .4  exp 

j(xsin<pR) 

=  Rex  p 

-  —  7'(*sin^r) 

_  vo 

v0 

V,n 

(119a) 


CO  ,  . 

0)  , 

CO 

Z0  cos  4 p,H  exp 

- j(xsm<Pj) 

-Z0  cos  cpRA  exp 

- j(xsm(pR) 

=  Zm  cosi/yRcxp 

- j(xsm(pT) 

L  vo 

L  vo 

L  Vm  \ 

(119b) 

The  foregoing  equations  are  identical  to  equations  and  hence  we  have  for  negative  refraction. 

Solution  for  reflection  in  a  negative  s  and  negative  ju  material  off  an  air  interface 

Based  on  what  we  learned  in  the  previous  cases  we  can  readily  write  down  the  equations  for 
figure  3. 


H  .  =H  exp  —  j(x  sin  (p,  +  zcos^,) 

V,n 

(120a) 

Ex  j  =  Zm  cos  (P[H  exp  —  j(x  sin  cp,  +  z  cos  (p, ) 

V,n 

(120b) 

Ez  j  =  —Zm  sin  cpjEl  exp  —  j(x  sin  (pt  +  z  cos  q>j ) 

v 

m 

(120c) 

H  R  =  Aexp  —  j(xsin(pR  - zcoscpR ) 

V 

m 

(121a) 

ExR=  — Z  m  cos  cpRA  exp  —  j(x  sin  <pR—z  cos  <pR ) 

v 

m 

(121b) 

EzR=  — Z  m  sin  <pRA  exp  —  j(x  sin  <pR—z  cos  cpR ) 

v 

m 

(121c) 

H  T  =  B  exp-  —  j  (-x  sin  (pT  +  z  cos  (pT  ) 
v0 

(122a) 

Ex  T  =  Z0  cos cpTB  exp-  —  j(-x sin <pT  +  z  cos <pr) 

vo 

(122b) 

E,t  =  Z0  sin^»rRexp-  —  j(-xsm<pT  +  z  cos  cpT) 

(122c) 

Matching  boundary  conditions  at  z  =  0  gives 
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(123a) 


H  exp  —  j(x sin <p,)  +  A  exp  —  j(x sin <pR)  =  B  exp  —  j(x sin cpT ) 

V,„  v„  vn 


Zm  cos cpjH  exp  —  j(x sin (p,)  -  Zm  cos (pRA  exp  —  j(x  sin (pR)  =  ZQ  cos cpTB  exp  —  j(x sin cpT ) 

Vm  Vm  V0 

(123b) 


The  solution  is 


<PR=<Pi  (124a) 

sin  <Pi  =  sin  <pT  (1 24b) 

H  +  A  =  B  (124c) 

Zm  cos  (p,(H  -  A)  =  Z 0  cos  (pTB  =  Z0  cos  (pT  {H  +  A)  ( 1 24d) 


A_  Zm  cos  <p,-Z0  cos  <pr  r 
Zm  cos  cp,  +  Z0  cos  cpT 

2z„,cos  <Pj_  H 
Zm  cos  (pI  +Z0cos^r 


(124e) 

(124f) 


4.  Conclusion 


We  derive  a  system  of  propagation  equations  and  model  in  a  Double  Negative  material  in  a  way 
that  differs  from  previous  derivations.  Our  derivation  is  based  entirely  of  on  the  idea  that  real 
materials  always  have  some  loss,  and  because  of  this,  wave  energy  traveling  in  a  certain 
direction  must  always  be  accompanied  by  a  loss  of  energy  in  that  direction.  Additional 
mathematics  is  not  required.  Energy  losses  per  unit  length  of  travel  are  finite,  and  can  be 
extremely  small.  Propagation  in  a  lossless  DN  media  is  found  as  the  mathematical  limiting 
solution  of  an  extremely  small  energy  loss  per  unit  length.  Our  system  of  equations  and 
methodology  is  used  to  derive  the  reflection  and  refraction  equations  between  a  positive  s(co)  and 
ju(co)  and  a  DN  material.  These  results  are  in  agreement  with  other  predictions. 
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